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\^ ', Abstract 

(N 

1^ . We study the sensitivity of our recent MSTW 2008 NLO and NNLO PDF analyses to 

^P , the values of the charm- and bottom-quark masses, and we provide additional public PDF 

sets for a wide range of these heavy-quark masses. We quantify the impact of varying 
rric and mi, on the cross sections for W, Z and Higgs production at the Tevatron and 
the LHC. We generate 3- and 4-flavour versions of the (5-flavour) MSTW 2008 PDFs by 



;h ' evolving the input PDFs and as determined from fits in the 5-fiavour scheme, including 

the eigenvector PDF sets necessary for calculation of PDF uncertainties. As an example 
of their use, we study the difference in the Z total cross sections at the Tevatron and LHC 
in the 4- and 5-flavour schemes. Significant differences are found, illustrating the need to 
resum large logarithms in Q"^ /m^ by using the 5-flavour scheme. The 4-flavour scheme 
is still necessary, however, if cuts are imposed on associated (massive) 6-quarks, as is the 
case for the experimental measurement of Zhh production and similar processes. 



1 Introduction 

Parton distribution functions (PDFs) determined from global analysis [IHS] of fixed-target and 
collider data, mainly from deep-inelastic scattering (DIS), are an essential input to theory 
calculations for hard processes at hadron colliders such as the Tevatron and LHC In addition 
to the fitted input PDFs, several other parameters enter into the global fits, which affect both 
the PDFs obtained and predictions for hadronic processes. One important example is the value 
of the strong coupling as and its uncertainty [OHH] • Others, which are the focus of the present 
paper, are the values of the heavy-quark masses and the scheme used to calculate heavy-quark 
contributions to observables. In particular, while the precise values taken for tjIc and irth may not 
be crucial for the processes included in the global PDF analyses, calculations of processes with 
explicit c- and 6-quarks might be more sensitive to these values, and it is therefore desirable to 
have PDFs available which have been consistently fitted and evolved with corresponding values 
of TTic and r/if). 

We will study two topics, both of which concern the treatment of the heavy charm and 
bottom quarks in global parton analyses. The first topic, which is the subject of Section |3l 
concerns the sensitivity of the MSTW 2008 global parton analysis [1] to the values of the heavy- 
quark masses m/j, with h = c,b. In Ref. [I] these masses were taken to be rric = 1.40 GeV 
and rrih = 4.75 GeV. Here we perform global fits for a range of rric and rrih about these values, 
with asiM^) allowed to be a free parameter. In this way, we determine the values of rric 
and TTib preferred by the data, and the correlations between these values of rrih and asiM"^). 
Due to a significant correlation between rric and as{M'^), we also perform fits with varying 
rric but with as^M"^) held fixed. We study how the cross sections for W, Z and Higgs boson 
production at the Tevatron and the LHC depend on the choice of rric and rrib, and we provide a 
recommendation on how to include the uncertainty arising from ruh in a general cross-section 
calculation. 

The second topic, described in Section HI is the generation of sets of 3- and 4-fiavour PDFs 
corresponding to the MSTW 2008 pQ 5-flavour sets of PDFs. We follow the same procedure pre- 
viously used to generate the 3- and 4-flavour versions [S] of the 5-flavour MRST 2004 PDFs |10j . 
i.e. the input PDFs (and as) aX Ql = 1 GeV^ determined from fits in the 5-flavour scheme are 
used as the initial condition for 3- or 4-flavour evolution. However, going beyond Ref. [9], we 
also provide 3- and 4-flavour eigenvector PDF sets to enable calculation of PDF uncertainties, 
and we provide 3- and 4-flavour PDF sets for a wide range of rUc and mi, values, respectively. 
As an example of the use of the central 4-flavour PDFs for the default quark masses, we com- 
pare the total cross sections for Z production at the Tevatron and LHC in the 4- and 5-flavour 
schemes. However, first we begin with a short resume of the alternative "schemes" to treat 
heavy flavour^j, and in particular, an explanation of what precisely we mean by 3-, 4-, and 
5-flavour PDFs. 



^Strictly speaking, it would be better to say alternative "techniques", since the use of the word "scheme" is 
usually reserved for an alternative choice in the ordering of the perturbative expansion, or a particular separation 
of contributions between the coefficient functions and parton densities — an ambiguity inherent in QCD. 



2 Schemes for the treatment of heavy flavours 

It is appropriate to briefly recall the various schemes for the treatment of heavy flavours in 
global parton analyses. In PDF analyses it is common to start evolving upwards from Q^ = 
Qo ~ 1 GeV^ with the distributions of the three light quarks {u,d,s), assuming that they 
are massless. As we evolve upwards, we have the choice to turn on the heavy-quark (c, b, t) 
distributions as we pass through their respective transition point, for which a convenient choice 
is fi'jp = m\. As we pass through a transition point, the number of active quarks is increased 
from n to n + 1, and the parton densities are related to each other perturbatively, i.e. 

k 

where the perturbative matrix elements Ajk{^'^p/in\) contain ln(/i|,/m^) terms which are known 
to C(a|) [H] and C(a|) [12]. The "x" arguments have been suppressed in Eq. ([1]), and the 
symbol ® is shorthand for the convolution 

/■^ Ax' 

f®g^ —f{x')g{x/x'). (2) 



X 

Eq. ([1]) relates the fj^ifJ^'p) set of partons to the f^^^ifJ^'p) set, guaranteeing the correct evolution 
for both the n and n + 1 regimes. We make the simplest choice, fi'j? = Q"^, for the factorisation 
scale. 

Hence, we have to decide whether or not to keep a heavy quark as just a final-state particle, 
and not as a parton within the proton. We may choose to keep just the three light flavours as 
parton distributions. We will call this the 3-flavour scheme (3FS), though it is often referred to 
as the fixed flavour number scheme (FFNS). Alternatively, we may include the c-quark in the 
evolution above Q^ = m^ and generate 4- flavour PDFs in a 4- flavour scheme (4FS). Actually, 
in the global MRST/MSTW parton analyses we also include the 6-quark distribution in the 
evolution above Q^ = ml, but not the t-quark above Q^ = mf, so we generate 5- flavour sets of 
parton distributions in a 5-flavour scheme (5FS). So to be precise, in our definition of n/-flavour 
parton sets, nf refers to the maximum number of quark flavours in the evolution. 

In each rij-flavour scheme (n^FS) the structure functions are given by the usual convolution 
of coefficient functions and parton distributions: 

F{x,Q') = J2cr''iQVrnl) ® f-{Q% (3) 

j 

where the sum j is over the gluon and the (variable, depending on Q"^) number of active quark 
flavours, n < Uf. We have a choice in how to choose Uf and define the coefficient functions. 
One simple choice is to fix n = Uf = 3. For the heavy flavours, all the rrih dependence then 
occurs in the coefficient functions, and these are called the FFNS coefficient functions. The 
structure functions may be written in the form 

F(x,Q2)= Yl Cf'\QVml)0ffiQ^). (4) 

j=u,d,s,g 



However, the sum of the a|ln"^(Q^/m|) terms, with m < k, is not included in the perturbative 
expansion. Thus the accuracy of the expansion becomes increasingly unreliable as Q^ increases 
above m\. In addition, there is the problem that the full mass dependence of the coefficient 
functions is known up to NLO |13], but is not completely defined at NNLO, i.e. the a| coeffi- 
cient, 6*2 f^' , for F2 is not fully known, see Ref. [12]. (Here, the outer subscript ''hg" denotes 
the g ^ h contribution to the heavy-fiavour structure function F2.) 

As an aside, we note that it would be possible to treat the charm quark as light, and the 
bottom quark as heavy. Then it would be possible to express the structure functions and 
cross sections in terms of four light quarks with all the mass dependence of the bottom quark 
contained in the coefficient functions. This is sometimes called the 4-flavour FFNS. However if 
one needs to use scales Q^ < ml and in the vicinity of ml, as for example in the global PDF 
analysis, then the charm-quark mass dependence, and the charm transition point (if Q^ < m^), 
should be included. Thus, this 4-flavour FFNS is only applicable in a restricted range of Q^, 
and it will not be considered further here. 

The alternative, and better, approach is to use the 4-flavour (and 5-flavour) PDFs of the 
variable flavour number schemes (4FS, 5FS). The simplest variable flavour number evolution 
procedure is to evolve treating some (or all) of the heavy quarks as massless, but to turn on 
the distributions at the appropriate transition points, Q^ = m,\. That is, to assume that these 
heavy-quark distributions evolve according to the splitting functions for massless quarks. Thus, 
the resummation of the large logarithms in Q"^ /m,\ is achieved by the introduction of heavy- 
flavour parton distributions and the solution of the evolution equations. It is motivated by 
the observation that at high scales the massive quarks behave like massless partons, and the 
coefficient functions are simply those in the massless limit, e.g. for structure functions 

F(x,Q^) = 5^c;^™^^®/7(g^), (5) 

j 

where, as in Eq. ([3]), the sum j is over the gluon and the (variable) number n < nj oi quarks 
that are active during the evolution. This is the so-called zero-mass variable flavour number 
scheme (ZM-VFNS). However, the ZM-VFNS has the failing that it simply ignores 0{m\/Q'^) 
corrections to the coefficient functions, and hence it is inaccurate in the region where Q^ is not 
so much greater than m\. Thus, strictly speaking, the ZM-VFNS does not define a scheme, 
but rather an approximation. 

So we have two approaches where the treatment of heavy flavours is relatively simple. 
The 3-flavour scheme (or FFNS), appropriate to the region Q^ < m^, and the ZM-VFNS, 
appropriate to the region Q"^ ^ m^. Clearly, for precision parton analyses, we must use a 
so-called general-mass variable flavour number scheme (GM-VFNS) which smoothly connects 
these two well-deflned regions, so as to reduce to the (3-flavour) FFNS in the low Q^ limit and 
to the ZM-VFNS in the high Q^ limit (up to possible terms of higher order in as)- In particular, 
in Eq. ([3]) the value of n increases by one each time we reach Q^ = ml in the evolution, and 
the coefficient functions C"*^^^^((5^/m^) interpolate smoothly from Cj ' (Q^/m^) to Cj ^ , 



for the maximuin value of n = nf. There is some freedom in how one does this, and in the 
MSTW 2008 analysis pp we use the definition of the GM-VFNS as described in Ref. [Il]@ We 
also note that actually a 5-fiavour GM-VFNS is used, since we do not include the top quark 
in the evolution. Note that this means that processes such as associated production of Higgs 
bosons with top quarks should be calculated using the 5-fiavour PDFs with the full 2—7-3 
subprocess matrix elements, i.e. gg, qq — )■ Hti. In a 6FS, one would calculate this via ti — )■ H, 
introducing a top-quark PDF. Although in principle this method would include resummed 
higher-order corrections ~ [a5ln(M|f/m^)]", we would expect that for all practical purposeqj 
at the LHC, the 5FS approach, gg, qq — )■ Hti, supplemented by NLO corrections, would be 
adequate, see also Section IT2l 

As just noted, the use of a GM-VFNS is nowadays essential for the determination of a set of 
precision PDFs from the datajj However, there are processes where the full mass dependence of 
the coefficient functions for c and b production processes are only known in the case where the 
heavy quark is treated as a final-state particle, and not as a parton in the proton. For example, 
FFNS parton distributions are needed for use with the HVQDis [19ii20j and mcOnlo [2T] 
programs. Thus in this paper we make available 3- and 4-fiavour PDFs. We also give an 
example of their use. This is the subject of Section HI 

3 Dependence of PDFs on the heavy-quark masses 

In the MSTW 2008 global parton analysis pLj we used rric = 1.40 GeV and m^ = 4.75 GeV, 
where these are the pole mass values. Some limited justification for choosing these values was 
briefiy given at the end of Section 3 of Ref. [1]. To summarise, there was little data constraint on 
nib, so it was simply fixed at a value of 4.75 GeV, close to the calculated MS mass transformed 
to the pole mass. The fixed value of rric = 1.40 GeV was close to the best-fit value at NLO if 
treated as a free parameter, but a little higher than the best-fit value at NNLO, and lower than 
the calculated MS mass transformed to the pole mass. We discuss limitations in constraints on 
the pole masses by transforming from the calculated MS masses below in Section 13.31 

Here, using exactly the same data sets, we repeat the global analysis for different fixed 
values of rric and mi,, but with as{M'^) left as a free parameter. The MSTW 2008 global 
fit pj used a wide variety of data from both fixed-target experiments and the HERA ep and 
Tevatron pp colliders. Neutral-current structure functions {F2 and Fl) were included from fixed- 
target lepton-nucleon scattering experiments (BCDMS [221123], NMC [21I2S], E665 [26] and 



^At NNLO this involves some modelling of the 0{ag) FFNS coefficient functions which are not fully calcu- 
lated. In the GM-VFNS we only need these in the low-Q^ regime and we approximate them using the known 
results of the small-x 15j and threshold limits [TB] (see Refs. pTlfTB] for recent refinements). 

3Note that for Mh < 1 TeV, as \n{M]j / 771"^) < 0.35. 

^In fact we only use the GM-VFNS for structure functions in DIS, whereas hadron collider cross sections 
are calculated in the ZM-VFNS under the assumption that mass effects are negligible at the typically large 
scales. (Low-mass fixed-target Drell-Yan production is also calculated in the ZM-VFNS, but here heavy flavours 
contribute <C 1% of the total, so inaccuracies induced by using the ZM-VFNS are in practice irrelevant.) 
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(b) NNLO: 
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Table 1: Fit quality and as{Ml) for different m^ values at (a) NLO and (b) NNLO. 



SLAC [271 - [2^ ). low-mass Drell-Yan cross sections from the E866/NuSea experiment 
and charged-current structure functions {F2 and xF^) and dimuon cross sections from neutrino- 
nucleon scattering experiments (CCFR/NuTeV [32l[33] and CHORUS |3l]). From the HERA 
experiments, HI and ZEUS, data were included on neutral- and charged-current reduced cross 
sections (crj^^ and cr^*") [2SH13], the charm structure function (^2^) [HHSU] . and inclusive jet 
production in DIS [5TH53|. From the Tevatron experiments, CDF and D0, Run H data were 
included on inclusive jet production [5HI55] . the lepton charge asymmetry from W decays [561157] 



and the Z rapidity distribution [581I59]. A more detailed description of the treatment of each 
of these data sets can be found in Ref . [1] . Note that more precise HI and ZEUS data on F| 
(and the beauty structure function, F^, not included originally) are now available, but we stick 
to fitting exactly the same data sets as in the MSTW 2008 analysis. 



3.1 Dependence on charm-quark mass rric 



The sensitivity to the charm mass of the data in the global PDF analysis, and the subset of -F| 
data, in terms of the goodness-of-fit measure x^, is shown in Table [T] at NLO and NNLO. The 
global best-fit values of rric, found by varying rric in steps of 0.01 GeV, are rUc = 1.45 GeV at 
NLO and rric = 1-26 GeV at NNLO. Recall that in Section 3 of Ref. [Tj we instead quoted best- 
fit values of rric = 1-39 GeV at NLO and rric = 1-27 GeV at NNLO. In fact, the more detailed 
analysis performed here, by scanning rric in steps of 0.01 GeV, revealed a rather flat minimum 
of the global x^ versus rric, with a slight double-minimum structure at NLO (the two minima 
being at 1.39 GeV and 1.45 GeV). Hence the best-fit values we quote here, of rric = 1-45 GeV 
at NLO and rric = 1-26 GeV at NNLO, differ slightly from the best-fit values given in Ref. [1]. 
Note from Table [H^a) that there is a clear correlation between rric and the value of 05 (M|) 
obtained from the NLO fit, and that from Table [T](b) this correlation is reduced at NNLO. 

In Fig. [1] we show the NLO and NNLO global x^ profiles as a function of rric- We see that at 
NLO there is a preference for rric = 1-45 GeV, with a la uncertainty of Iq 06 GeV. The global 
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Figure 1: Dependence on rric of the global x^, relative to the value at the best-fit rric, at (a) NLO 
and (b) NNLO, fitting to the same data sets used in the MSTW 2008 global analysis [T]. The 
values quoted for rric correspond to treating as{M'^) as a free parameter, and the corresponding 
global x^ values are shown by the closed points (•) joined by the continuous curves. The 68% (la) 
and 90% confidence-level (C.L.) uncertainties on rric are indicated by the vertical dashed lines. The 
open points (o), joined by the dashed curves, are the global x^ values with a5(M|) held fixed. 



data prefer a lower value of rric at NNLO, namely rric = 1.26 GeV, with a la uncertainty of 
lo'o4 GeV. (For the higher values of rric, the NNLO fit considerably undershoots the moderate 
Q^ deep-inelastic data.) Here, the la uncertainty is obtained in exactly the same way as 
described for the uncertainty on as{M^) in Ref. j6j, i.e. by examining the x^ profile of each 
data set included in the global fit. The distinguishing power of the various data sets is shown 
in Fig. [2J Here, the points (•) indicate the values of rric for which the x^ fo^^ each data set is 
minimised (within the context of the global fit), while the inner error bars extend across the 
68% confidence- level (CL.) region and the outer error bars extend across the 90% CL. region 
(see Section 4 of Ref. [6] for the precise definitions of the 68% and 90% CL. limits). As may be 
expected, the -F|(a;, Q^) data of the HI and ZEUS collaborations are the most discriminating. 
Indeed, they dominate the determination of rric, as can be seen from Fig. [3], which shows x^ 
versus rUc for the F2 data alone. For the HERA inclusive data, changes in the value of ttIc are 
partially compensated by changes in the gluon and in as- The NMC data prefer a lower value 
of rUc giving a quicker evolution near threshold. Similarly the BCDMS data mainly prefer a 
lower rric, but only due to the correlation with a lower value of as- Supplementary plots of 
the x^ profiles versus rric for all data sets in the global fit are available from Ref. ^U\. The 
experimental uncertainty on rric, given in Fig. [T] and indicated by the horizontal dashed lines in 
Fig. [21 is chosen to ensure that all data sets are described within their 68% or 90% CL. limits. 

The dependence of the PDFs on the value taken for rUc is shown in Fig. |H The plots 
show the ratio of the NLO PDFs obtained with rUc = 1.3 GeV and rric = 1.5 GeV to those 
of the MSTW 2008 analysis (which assumed rric = 1-4 GeV) at scales of Q^ = 4 GeV^ and 
Q2 = 10^ GeV^ We see that the ratios lie well within the 90% CL. PDF uncertainty. The 
exception is the charm distribution at Q^ = 4 GeV^ where, as expected, the different charm 
transition points have an appreciable effect. 

3.2 Dependence on bottom-quark mass rrib 

We repeated the whole exercise for rrib- That is, we performed a series of global fits for different 
values of mi,- The results are shown in Table [2J Now there is less sensitivity to the value of 
the quark mass, and essentially no correlation between mj, and a5(M|). Indeed, the global 
X^ stays fairly flat all the way down to nib = 3 GeV. For the lower values of rrib, there is a 
slightly better description of the HERA data, including F2{x,Q'^)- A similar conclusion holds 
at NNLO, but with about half the change in global x^j that is, even less sensitivity to rrif,- 

Note that the HERA F|(x, Q^) data [I51li71lgTtiB5] are not included in the global fits, neither 
in MSTW 2008 [Ij, nor in the fits described above. In Fig. [5] we compare the predictions of the 
NNLO fits with varying rrib values to the F2 data from HI [61]. We also show a few of the (less 
precise) ZEUS data points |62]. In Table [2] we give the x^ for the HI F2 data accounting^ for 
all 24 sources of correlated systematic uncertainty, and we also give the simple addition with 



^We use Eqs. (38-40) of Ref. [T], noting that there is a typo in Eq. (40) of that paper and cr""''"'''' should 
appear squared in the expression for Akk' ■ 
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Figure 2: Ranges of rric for which data sets are described within their 90% C.L. limit {outer error 
bars) or 68% C.L. limit {inner error bars) in the (a) NLO and (b) NNLO global fits, with a5(M|) 
as a free parameter. The points {•) indicate the values of rric favoured by each individual data set, 
that is, the values for which the x^ is minimised. The experimental uncertainty on rric, indicated by 
the tiorizontal dashed lines, is chosen to ensure that all data sets are described within their 68% or 
90% C.L. limits defined by Eq. (2) of Ref. [6J. 
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Figure 3: Dependence on nic of the x foi" HERA data [44l - [50] on the charm structure function, F2, 
relative to the value at the global best-fit nic, at (a) NLO and (b) NNLO. The dosed points (•) have 
a5(M|) as a free parameter, and are fitted to a quadratic function of mc shown by the continuous 
curves, while the open points (o, mostly hidden) have as{M'^) held fixed at the MSTW 2008 value. 
The horizontal dastied lines in the plots indicate the 68% and 90% C.L. limits, determined according 
to a "hypothesis-testing" criterion, see Section 4 of Ref. [6j. 
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Table 2: Fit quality and as{M^) for different m,, values at NLO and NNLO. Note that the F| 
data [61] are not included in the global fit. 
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Figure 4: rric dependence of gluon, singlet and charm distributions at NLO at two different Q^ 
values, 4 GeV^ (left) and 10'^ GeV^ {right), compared to the 90% C.L. PDF uncertainties. 
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Figure 5: Comparison of the HERA F2^(x, Q^) data [SlEl] with NNLO predictions obtained using 
PDFs from global fits with different values of nih. Note that the F^ data are not included in these 
global fits. The curves, in decreasing order, correspond to nib = 4, 4.5, 4.75, 5, 5.5 GeV. 
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Figure 6: nib dependence of bottom-quark distribution at NLO at two different Q^ values, 
(a) 40 GeV^ and (b) 10^ GeV^, compared to the 90% C.L. PDF uncertainties. 

the global x^- We see that the F2 data show a slight preference for nib ~ 4.75 GeV at NLO and 
rub ~ 5.00 GeV at NNLO. The global fits (including the F^ data, obtained by simply adding 
the x^ values), would prefer rub ^ 4.50-4.75 GeV at NLO and nib ~ 4.50-5.00 GeV at NNLO, 
although admittedly the x^ profiles are quite fiat, see Table [21 We conclude that the global 
data do not meaningfully constrain nib, and that, in view of Table [2] and Fig. [5], the choice 
rrib = 4.75 GeV made in the MSTW 2008 [1] analysis is completely satisfactory. 

The dependence of the bottom-quark distribution on the value taken for nib is shown in 
Fig. [6l The plots show the ratio of the NLO PDFs obtained with rUb = 4.25 GeV and rUb = 
5.25 GeV to those of the MSTW 2008 analysis (which assumed nib = 4.75 GeV) at scales of Q^ = 
40 GeV^ and Q^ = W GeV^. We see that the 6-quark ratios lie well outside the 90% C.L. PDF 
uncertainty at Q^ = 40 GeV^ and even at Q^ = 10^ GeV^. The gluon distribution and the 
other quark flavours are hardly affected by the value of nib (not shown here). 



3.3 Pole masses of heavy quarks 

Note that in our analyses m^ refers to the pole mass of the charm quark. Both the transition 
matrix elements defining the boundary conditions for heavy-quark evolution [11] and the heavy- 
quark coefficient functions used to define the GM-VFNS [TT] are defined using the on-shell 
renormalisation scheme for the heavy quark. Unlike asin'^), the pole mass rric is, in principle, 
a physically-defined (spoilt by confinement which gives power corrections, see below) quantity 
which is stable to the order of perturbation theory. So the values that we obtain in the NLO 
and NNLO analyses should, in principle, be the same. However, we see from Fig. [1] that the 
values are only marginally compatible. This is due to a variety of reasons. The perturbative 
expansion of heavy-flavour coefficient functions is not as convergent as one might hope, with 
large corrections near threshold and in the small-x limit. So even though the pole mass might, 
in principle, be the same quantity at two different perturbative orders, the calculation at lower 
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orders is missing significant corrections at higher orders which will affect the value of the pole 
mass obtained in an extraction from data. Perhaps most important in this respect is the large 
negative boundary condition (at Q^ = m|) for the heavy-flavour distribution at NNLO, which 
has the effect of flattening the slope dF2/dlnQ'^ (see Ref. [H] for details) and implies a lower 
mass to correct for this. There is also more potential variation in the GM-VFNS definition 
at NLO than NNLO, and there would be slightly more stability in preferred pole mass in, 
for example, the "optimal" scheme introduced in Ref. [Hlj (although this does not remove the 
variation and is the subject for a future study). Finally, as briefly noted above and discussed 
further below, the definition of the pole mass is contaminated by non-perturbative corrections. 
These are reflected in the divergence of the perturbation theory and power corrections are, in 
practice, difficult to disentangle from higher-order corrections, so extractions at different orders 
will have a variation from this source. We now examine this issue in more detail. 

Let us consider the values we extract and use as default. In the Review of Particle 
Physics [65], the authors quote values for the ^^running" masses, rricifi = rric) and mh{mb), 
in the MS scheme: 

m,{fi = m,) = 1.27 +°:?^ GeV, m,{m,) = 4.20 +°:J^ GeV, (6) 

which they evaluate from a whole series of independent determinations. This implies a precise 
determination of massesjj The best precision of the pole mass should be obtained by using this 
MS value (the best determinations coming from NNNLO calculations) with the most accurate 
conversion factor, also at NNNLO. However, the conversion of the ^^running" mass to the pole 
mass is problematic, since it relies on a very weakly convergent perturbative expansion. Using 
our NNLO value of as, the term- by-term conversion factor (obtained from the formula in 
Ref. [65]), starting at zeroth order and going to the highest known order of a| [681169] . for the 
bottom-quark mass is (1 + 0.095 + 0.045 + 0.036) and for charm is (1 + 0.16 + 0.14 + 0.18)0 The 
former implies that convergence is ceasing, while the latter has no real convergence at all, since 
it relies on a much higher value of as (and larger coefficients due to less cancellation between 
individual terms as the number of light quark flavours decreases [681169] ). 

Indeed, the contamination by renormalons of the perturbative series for the transformation 
of a running quark mass in the MS scheme to the pole mass has been known for a long time [7H 
[72]. These papers estimate the intrinsic ambiguity of the perturbative series due to infrared 
contributions, translating into an uncertainty on the pole mass of at least 0.05 GeV, but more 
likely of order 0.1-0.2 GeV. The best estimate one can obtain from the perturbative series 



^There are numerous individual determinations of both charm- and bottom-quark masses in the MS scheme 
which are actuahy even more precise than the above — for recent examples see Refs. [661E3- However, these 
can vary from each other by considerably more than the quoted uncertainties, so in the Review of Particle 
Physics [65j account is taken of the spread, as well as the individual uncertainties. 

^Use of an NNNLO value of as would make little difference. As shown in Ref. [70 , NNNLO extractions 
of as{M'^) are very similar to the corresponding NNLO values, and the slightly quicker running at NNNLO 
would marginally increase as at the scale of the charm- and bottom-quark masses, i.e. the convergence of the 
perturbative series would likely be very slightly worse. 
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is usually defined by including all terms in the series until they cease to fall in magnitude 
from one order to the next. For the bottom quark, under the plausible assumption that the 
unknown 0{ag) term is similar in size to the calculated 0(«|) term, this would give a pole 
mass of rub = 4.9 GeV. The uncertainty on this value from the conversion is either the estimate 
from Ref. [72] or roughly the size of the last term included in the series, i.e. 0.15 GeV. The 
good correspondence of the two estimates of the uncertainty is expected, and suggests that the 
series is indeed truncated at the correct point. If this uncertainty is added (approximately) in 
quadrature with the uncertainty on the MS scheme mass, Eq. (Q, we obtain mf, = 4.9±0.2 GeV. 

The determination of the charm-quark pole mass by conversion from the MS scheme mass 
is more problematic since, as noted above, the perturbative series displays no clear convergence 
at all. However, as discussed in Ref. [75], the leading renormalon contribution is the same for 
different masses, i.e. Sm'^ oc Aq^j-,, independent of quark fiavour. This means the difference 
between rub and ttIc can be well determined, and turns out to be 3.4 GeV with a very small 
uncertainty [73]. Using this result we obtain a best determination of the charm-quark pole 
mass of TTic = 1.5 ± 0.2 GeV. There is an uncertainty on rrib — rric due to the fact that the non- 
dominant renormalon contributions, 5m? oc Aq^j-j/m^ and beyond, do not cancel. However, the 
uncertainty on rUb — rric is small compared to the 0.2 GeV uncertainty on rub and hence does 
not affect the uncertainty on rric at the quoted accuracy when added in quadrature. 

These pole mass values obtained from MS conversion: 

rric = 1.5 ±0.2 GeV, nib = 4.9 ±0.2 GeV, (7) 

are slightly high compared with our default values (1.4 GeV and 4.75 GeV) and our best-fit 
values of rric, particularly at NNLO, but even in the most discrepant case the error bars overlap. 
Combining our best-fit rric values and the pole mass determination given in Eq. (ITj), it seems that 
a range for rric of about 1.25-1.55 GeV at 68% C.L. and 1.15-1.65 GeV at 90% C.L. would seem 
reasonable at present. With improved HERA data on F2 to come [71] this range can hopefully 
be narrowed in the near future. Turning now to the bottom mass, the limited information 
that we obtain from comparison to data, and from the determination of the pole mass given in 
Eq. (ITj), suggest that a range of m;, from about 4.65-5.05 GeV at 68% C.L. would seem sensible. 
As we will see in the next section, LHC cross sections that do not involve explicit 6-quarks are 
not very sensitive to the PDF variation with rrib, so using our grids at 4.5 and 5 GeV should give 
a good estimate of the uncertainty due to rrib at 68% C.L., and the grids at 4.25 and 5.25 GeV 
should give a conservative estimate of this source of uncertainty at 90% C.L. As with charm, 
improved HERA data on F2 are likely to limit the allowed spread in future. 

If we reach the stage where we become confident that the data are starting to constrain 
the pole masses to an accuracy clearly better than the renormalon ambiguity, it may become 
preferable to transform to a scheme where MS definitions are used instead, even though the 
mass is less directly related to a physical variable in this case. This would be appropriate 
if we become strongly constrained by data with Q^ ~ ml and data close to the kinematic 
threshold, W"^ = Q^iXjx — 1) = 4m^ (for neutral-current DIS), where non-perturbative effects 
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and the interplay between the leading-twist power-series and the power corrections become very 
important, and the ambiguities may be reduced in a different renormalisation scheme. Indeed, 
we have already noted in Section 9.2 of Ref. pQ that the lowest-Q^ EMC data on F^ [75] imply 
a non-perturbative correction to the cross section. 

3.4 Impact on VF, Z and Higgs production at the Tevatron and LHC 

Tables [3] and H] show how the predictions for the "standard candle" NNLO W and Z cross 
sections at the Tevatron and LHC change using the PDFs obtained when the charm- and 
bottom-quark masses are varied in the global fit. The cross sections are calculated in the 5- 
flavour ZM-VFNS as described in Section 15 of Ref. pTj, e.g. using the PDG 2008 |65] electroweak 
parameters. These changes are not due primarily to the changes in the heavy-quark PDFs 
themselves, which would in any case be in the opposite direction. Rather they are due to 
the changes in the light quarks, which evolve slightly more rapidly at small x to compensate 
for the slower turn-on of the heavy-flavour contribution to the structure functions when the 
heavy-quark masses are increased. In the free coupling case this is achieved mainly by an 
increase in as-, while for fixed coupling it occurs from an increase in the small- a; input gluon 
distribution. Additionally, the input sea quarks, mainly the less well-constrained strange quark 
distributions, also increase with increasing irth to compensate for the suppressed heavy-flavour 
contribution to small-x structure functions. At the Tevatron, varying the charm mass by ±0.15 
GeV from its default MSTW 2008 value leads to fairly small (0(1%) or less) changes in the W 
and Z cross sections, with heavier values of the charm mass giving slightly larger cross sections. 
The effects are a little more pronounced at the LHC: the cross sections vary by about ±2% 
for a/s = 14 TeV. These changes reflect the behaviour of the up-quark distribution shown in 
Fig. El^a) at the relevant x values (indicated). The sensitivity to the value of m^ of the W and 
Z cross sections at the LHC has previously been noticed in Ref. [76] . 

In Tables [3] and m we also give the changes to the cross section for the production via gluon- 
gluon fusion through a top-quark loop of a Standard Model Higgs boson of mass 120 GeV, 
calculated in the 5-flavour ZM-VFNS as described in Section 6.2 of Ref. [6]. From Table [3] we 
see that the variation of the Higgs cross sections with irtc is different depending on whether 
as is (a) allowed to vary or (b) kept fixed. This can be understood from the behaviour of the 
gluon distribution as rric is varied, see Fig. [7](b). When as is fixed, the change in the Higgs 
cross section reflects the change in the gluon distribution at the relevant x values (indicated). 
Thus the {14 TeV LHC, 7 TeV LHC, Tevatron} Higgs cross sections are {correlated, uncorre- 
lated, anticorrelated} with tti^. On the other hand, for variable as-, the additional correlation 
of as with rric (see Tabled]) results in {14 TeV LHC, 7 TeV LHC, Tevatron} Higgs cross sec- 
tions that are {correlated, correlated, almost uncorrelated} with ?7ic, i.e. at the Tevatron the 
anticorrelation of the gluon distribution largely cancels the correlation of as- 

From Table Hj we see that the W -, Z and Higgs production cross sections are much less 
dependent on the value of mj,. As seen from Table E] the correlation between mj, and a5(M|) 
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Table 3: (a) Difference (in percent) between the predictions for the W, Z and Standard Model 
Higgs {Mh = 120 GeV) NNLO production cross sections at the Tevatron and LHC, calculated using 
the variable-mc PDF sets, and the standard MSTW 2008 NNLO predictions with rric = 1.40 GeV. 
(b) Also shown are the corresponding results for variable rric with fixed as{M'^). 
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Figure 7: rric dependence of the NNLO (a) up-quark and (b) gluon distributions at the scales 
relevant for Z and Higgs boson {Mh = 120 GeV) production, respectively. Also shown is the rUc 
dependence when as{M'^) is held fixed. The values of a; = Mz,h/\/s probed at central rapidity 
(and p^'^ = 0) are indicated for the Tevatron {^/s = 1.96 TeV) and LHC (7 TeV and 14 TeV). 
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+0.1 

+0.2 

+0.3 



Table 4: Difference (in percent) between the predictions for the W , Z and Standard Model Higgs 
{Mh = 120 GeV) NNLO production cross sections at the Tevatron and LHC, calculated using the 
variable-mfc PDF sets, and the standard MSTW 2008 NNLO predictions with m,, = 4.75 GeV. Here 
as{Mz) is a free parameter in each fit, although the correlation with mf, is negligible, see Tabled 



is negligible, so we do not explicitly show results with fixed 05 (M|) since they would be almost 
the same as in Table HI In all cases, for both varying-rrtc and varying-rrift, the cross-section 
ratios cr^/a^ (and also cr'^ /a^ at the LHC) are almost completely unaffected. 

In order to assess the total uncertainty on a hadron collider cross section arising from varia- 
tions in PDFs, as and the heavy-quark masses, we need to define a prescription for combining 
the latter with the former. Our recommendation is to vary nic in the range rric = 1.40±0.15 GeV 
at 68% C.L. and ±0.25 GeV at 90% C.L., centred on the default value for fixed as, and to vary 
TUb in the range m^ = 4.75 ± 0.25 GeV at 68% C.L. and ±0.50 GeV at 90% C.L., and then to 
separately add the resulting cross-section variations in quadrature with the usual "PDF±a5'" 
uncertainty. Of course, this prescription does not account for all possible correlations between 
PDFs, as and mc,b, but it should be a sufficiently good approximation. The range of rric values 
is based on the slightly contrasting pulls of the pole mass determination from the MS conversion 
and our fit. As an example, in Table [5] we show the increase in uncertainty on the W, Z and 
Higgs total cross sections due to rric and mf, variations, compared to the "PDF only" [1] and 
"PDF±a5" [6J uncertainties. 

In Table [6] we compare the values of the heavy-quark masses, rric and rrib, used in the 
MRST/MSTW analyses, with the values taken in the most recent public PDF fits of other 
groups. The other two global PDF fitting groups, CTEQ [2111] and NNPDF 0, do not attempt 
to quantify the uncertainty coming from rric^, neither do the fits of "dynamical" PDFs by 
GJR/JR [77l[78]. The HERAPDFl.O fit [79] uses the same central values as MSTW08 [1], but 
additional PDF fits are also provided with rric = {1-35,1.65} GeV and rrih = {4.3,5.0} GeV, 
and this additional model uncertainty is recommended to be added in quadrature with the other 
uncertainties. The ABKM09 analysis [80] uses fixed values of rric = 1-5 GeV and rrib = 4.5 GeV 
to determine the central fit. However, additional pseudo-measurements oi rricfi are then added, 
with values given in the last line of Table O and rric and rrif, are taken as free parameters to 
calculate the covariance matrix used for the final error propagation. This means that each of 
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Tevatron, ^/s = 1.96 TeV 




Bi+i- ■ (7 


a^ 


Central value 


2.747 nb 


0.2507 nb 


0.9550 pb 


PDF only uncertainty 
PDF+q;5 uncertainty 


+ 1.8% 
-1.5% 
+2.2% 
-1.7% 


+1.9% 
-1.6% 
+2.2% 
-1.8% 


+3.1% 
-3.3% 
+5.4% 
-4.8% 


PDF+a;5+r7ic,b uncertainty 


+2.3% 
-1.8% 


+2.3% 
-2.0% 


+5.6% 
-5.1% 




LHC, v^ = 7 TeV 


R ^W 


Be+e- ' 0" 


a^ 


Central value 


10.47 nb 


0.958 nb 


15.50 pb 


PDF only uncertainty 
PDF+a^ uncertainty 


+ 1.7% 
-1.6% 
+2.5% 
-1.9% 


+1.7% 
-1.5% 
+2.5% 
-1.9% 


+1.1% 

-1.6% 
+3.7% 
-2.9% 


PDF+as+rric^i, uncertainty 


+2.7% 
-2.2% 


+2.9% 
-2.4% 


+3.7% 
-2.9% 




LHC, yi = 14 TeV 


Be, ■ a^ 


Be+e- ■ a^ 


a^ 


Central value 


21.72 nb 


2.051 nb 


50.51 pb 


PDF only uncertainty 
PDF+a^ uncertainty 


+1.7% 
-1.7% 
+2.6% 
-2.2% 


+1.7% 
-1.6% 
+2.6% 
-2.1% 


+1.0% 
-1.6% 
+3.6% 
-2.7% 


PDF+as+mc,b uncertainty 


+3.0% 
-2.7% 


+3.1% 
-2.8% 


+3.7% 
-2.8% 



Table 5: NNLO predictions for W, Z and Higgs {M^ = 120 GeV) total cross sections at the Teva- 
tron, 7 TeV LHC and 14 TeV LHC, with PDF uncertainties only [1], with the combined "PDf+as" 
uncertainty [6j, then finally also including the uncertainty due to rric and rub. The 68% C.L. uncer- 
tainties are given in all cases. We take jj^r = fj^p = My/^z,h- 



PDF set 


m, (GeV) 


mb (GeV) 


MSTW08 [H this work] 


1.40 ±0.15 


4.75 ± 0.25 


MRST06 [2] 




1.43 


4.30 


CTIO [3] 




1.30 


4.75 


CTEQ6.6 m 




1.30 


4.50 


NNPDF2.0 [5] 




v^ 


4.30 


GJR08/JR09 [771[78 


1.30 


4.20 


HERAPDFl.O 


m 


1 40+°-^° 


A 7C-+0.25 
^•'C>-0.45 


ABKM09 [8Cr 




1.50 ±0.10 


4.50 ±0.50 



Table 6: Values of ttIc and nib used in various PDF fits. The la uncertainties are given only for the 
three PDF groups who attempt to account for uncertainties on nicb in their analyses. 



20 



the public eigenvector PDF sets will be associated with different values of rric and m^, but these 
values are not readily accessible. 



4 3-flavour and 4-flavour scheme parton distributions 

As well as looking at the variation in the parton distributions as a function of the heavy-quark 
masses, it is also interesting to consider the PDF sets obtained in the framework of a different 
maximum number of active quark flavours. Hence, in this section we will consider our PDFs 
when charm becomes an active parton but bottom does not — the 4-fiavour scheme — and when 
both the charm and bottom quarks only appear in the final state — the 3-flavour scheme. We 
have argued on various occasions (e.g. in Section 4 of Ref. [Ij) that the use of a GM-VFNS, 
i.e. with up to five active quarks (or even six if we include top), is preferable. However, there 
are cases where the cross section has only been calculated with finite mass dependence for the 
fixed flavour number case; see, for example, the HQVDis program [1^1120] for details of final 
states in heavy-quark production in DIS. For this reason we make available sets of 3FS and 
4FS distributions with a variety of both charm and bottom masses. 

4.1 Obtaining 3-flavour and 4-flavour scheme parton distributions 

It might be thought preferable to obtain these lower active quark number PDF sets from a fit 
performed using FFNS scheme coefficient functions. However, as argued in Refs. [HEI], it is 
not actually so obvious that this is the case. This is largely because rather few of the data 
sets included in a truly global PDF fit can be kept, even at NLO, in a fit using the FFNS, due 
to lack of the full coefficient functions (even charged-current DIS coefficients are not known 
with full mass dependence at order q;| except for Q"^ ^ m\ [S5j). Hence, the central values of 
the PDFs are likely to be influenced as much by the lack of data as by the change of scheme. 
However, it is also a consideration that the lack of resummation of the large logarithms in 
Q"^ /m\ potentially affects the stability of the fit compared to the presumably more stable GM- 
VFNS fit. Ultimately a correct GM-VFNS will provide results very similar to the FFNS near 
to the transition points Q^ = ml anyway, so we deem it best to simply obtain the 3FS and 4FS 
PDFs from the inputs for the full fits performed using the GM-VFNS. 

When obtaining the 3FS and 4FS PDFs it is vital to make a self-consistent definition of 
the strong coupling constant as- It is generally the case that coefficient functions calculated in 
the 3- or 4-fiavour schemes are made in a renormalisation scheme where the contribution of the 
heavy quark decouples and the coupling itself does not include the heavy quark as an active 
fiavour. On this basis we define the PDFs using this definition of the coupling in the splitting 
functions. It is certainly possible to use a different definition of the coupling, but this must be 
applied universally, in both the PDFs and the coefficient functions. As illustrated in Ref. [9], 
the error, made from not doing so, can be a few percent. Indeed, the change in the coupling 
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Q^ (GeV^) 



Figure 8: The NLO strong coupling as(Q^) versus Q'^ , in the default 5FS with the 90% C.L. ex- 
perimental uncertainty, and in the 4FS and 3FS, taking the same input value of ols{Q\ = 1 GeV^). 

obtained by altering the number of heavy flavours is quite dramatic. In our default NLO fit we 
start with a value asiQl) ~ 0.49128 (with Q^ = 1 GeV^), which for a variable flavour number 
results in a^^(M|) = 0.12018. Restricting ourselves to a maximum of four active flavours, the 
same boundary condition at Qq gives a^ (M'^) = 0.11490, and for three active flavours gives 
: 0.10809. This is illustrated in Fig. [HI and more clearly in Fig. M^a), which shows 



(3) 



the ratio of the 3- and 4-flavour as to the 5-flavour one. 

We make available both the 3FS and 4FS NLO PDF sets for the full variety of fits with 
varying charm and bottom masses discussed in the previous section, including the case of 
varying rric with fixed as{M'^). We do not vary the quark mass at LO since the quality of the 
fit is already poor at LO. However, we do provide the 3FS and 4FS PDFs for the default masses 
at LO since these may be useful in some Monte Carlo event generators. We also make available 
the full range of 3FS and 4FS sets at NNLO, although the application of these is distinctly 
limited at present due to the lack of cross sections with full mass dependence calculated at this 
order. We will comment more on this in the following section. The ratios of the 3FS and 4FS 
NLO PDFs to those in the 5FS are shown in Fig. [H at Q^ = iq^ QeV^. In both cases, the 
gluon is larger in the 4FS and, particularly, in the 3FS, due to the fact that it splits into a 
smaller number of quark-antiquark pairs; see also Fig. [9](b). The effect of the increased growth 
of the gluon distribution is countered exactly for the leading term in asiQ"^) ln(Q^/m^) by the 
quicker decrease of the coupling. This is illustrated in Fig. [9l where the change in the gluon 
distribution comparing different flavour numbers is, to a good approximation, the inverse of 
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Figure 9: Ratio to the 5FS values of (a) strong coupling as and (b) gluon distribution at x = 0.01, 
versus Q^ at NLO in the 4FS and 3FS, compared to the 90% C.L. experimental uncertainty bands. 
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Figure 10: Ratio of 3FS and 4FS NLO PDFs to the standard 5FS NLO PDFs, at Q^ = 10^ GeV^ 
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that for as- This behaviour is only violated by higher-order corrections, so the gluon-driven 
change in the light quarks at small x is rather minimal. At higher x, the onset of the smaller 
coupling in the lower flavour number schemes means a slower evolution and consequently larger 
light-quark distributions, as seen in Fig. [TOl 

Finally, we also make available new eigenvector PDF sets for the 3- and 4-fiavour PDFs 
using the default quark masses at LO, NLO and NNLO. These are evolved from the saved PDF 
parameters at Qq = 1 GeV^ for the default MSTW fits, i.e. the "dynamical" tolerance values, 
T = (AXgiQjjj^j)^/^, are those determined from the 5-fiavour fit pQ. For the first time, this will 
allow PDF uncertainties to be consistently included in a 3FS or 4FS calculation. We have not yet 
generated the additional eigenvector PDF sets with varying as values needed for the "PDF+as" 
uncertainty calculation [6], although these could be provided at a future point if necessary. 
Note, however, that 3FS or 4FS calculations generally have large higher-order corrections, and 
an associated large (renormalisation and factorisation) scale dependence, which is common 
in processes with multiple hard scales. The theory uncertainty due to neglected higher-order 
corrections is therefore likely to dominate over any PDF or "PDF+as" uncertainty. Precise 
calculations of the latter quantities are therefore relatively less important than in typical 5FS 
calculations where higher-order corrections are more readily available. 

4.2 Comparison of Z total cross-section predictions in 4FS and 5FS 

We have already noted in Section 2 that while the parton distributions of a nFS, with appro- 
priate coefficient functions, can give an adequate description of the production of an (n + l)th 
parton near threshold, the accuracy of the perturbative expansion becomes increasingly unre- 
liable as the scale of the hard scattering process, Q^, increases above m\. In this section we 
illustrate this by considering the total Z cross section at hadron colliders. In particular, we 
compare the cross-section predictions in the 5FS (= MSTW 2008 [1]) and 4FS defined pre- 
viously. Since in this case Q^ = M| ^ ml, we would expect to see quantitative differences 
between the predictions due to the resummed [as'log(M|/m^)]"' terms which are implicit in 

(5) 

the 5FS, via the evolution of the 6-quark PDF and a^ , but absent from the 4FS. For reasons 
which will become apparent below, we consider the Z cross section at both NLO and NNLO. 

In the 4FS the b quarks are not considered as partons in the initial state, but contribute to 
the Z cross section via real and virtual contributions which first appear at 0(a|) in perturbation 
theory. Therefore the only difference in the predicted cross sections at NLO are (i) explicit b- 
parton contributions b-\-b — t- Z{+g) and g + b{b) — )■ Z + b{b) which only contribute in the nj = 5 
scheme (5FS), (ii) small differences in the light-quark and gluon distributions arising from the 
slightly different evolution in the two schemes, see Fig. [TUl and (iii) the difference in the values 
of a'-s\M^) and af{M^), see Figs. [8]and|9];a), which affect the size of the NLO K-iactor. We 
would expect the Tevatron Z cross sections to be more similar than those at the LHC, since the 
6-quark contributions to the 5FS cross section are much smaller at the lower collider energy. 
We do not explicitly include PDF and/or as uncertainties on the calculated cross sections, as 
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Tevatron, ^/s = 1.96 TeV 


i?-<o(4FS)(nb) 


i?-4Lo(5FS)(nb) 


i?-4Lo(5FS,6)(nb) 


-0^ 


0.1989 
0.0413 


0.1990 
0.0436 


0.0012 
-0.0002 


total 


0.2402 


0.2426 


0.0010 




LHC, v^ = 7 TeV 


i?-<o(4FS)(nb) 


i?-o-|Lo(5FS)(nb) 


i?-a|Lo(5FS,6) (nb) 




0.7846 
0.1206 


0.8023 
0.1285 


0.0205 
-0.0020 


total 


0.9052 


0.9308 


0.0185 




LHC, v^ = 14 TeV 


S-<o(4FS)(nb) 


5.a^Lo(5FS)(nb) 


S-a|Lo(5FS,6)(nb) 


^0^ 


1.6922 
0.2303 


1.7545 
0.2465 


0.0656 
-0.0050 


total 


1.9225 


2.0009 


0.0601 



Table 7: NLO predictions for the total Z cross section (multiplied by leptonic branching ratio B) 
at the Tevatron and LHC using MSTW 2008 NLO PDFs [Ij as input, broken down into the ag 
{n = 0,1) contributions, with {q = u,d,s,c] a^^^; 4-flavour MSTW 2008 NLO PDFs} in the 4FS 
calculation and {q = u,d,s,c,b] aj^; 5-flavour MSTW 2008 NLO PDFs} in the 5FS calculation. 
The final column gives the contribution to the 5FS cross sections from processes where the Z couples 
directly to b quarks. 



these are presumed to be more or less the same for the 4FS and 5FS calculations. Contributions 
from top quarks are also not included, as these have been shown to be very small in Ref. [S3], 
and in any case should be added to both the 4FS and 5FS cross sections. 

The results at NLO are shown in Table [7]cl The 4FS cross section is smaller by 1.0%, 
2.8% and 3.9% at the Tevatron, LHC(7 TeV) and LHC(14 TeV) respectively. The C(a^) 
contributions differ by more than the 0{a^^) contributions, due to the ~ 5% differences in the 
asiM"^) values in the two schemes. Also shown in Table [7] (final column) are the contributions 
to the 5FS cross sections from the b + b ^ ^i+d) ^^^ 9 + b(b) -^ Z + b(b) processes, which 
evidently account for 0.4%, 2.0% and 3.0% of the total at the Tevatron, LHC(7 TeV) and 
LHC (14 TeV) respectively. Thus these contributions account for the bulk of the differences in 
the 4FS and 5FS cross sections at the LHC energies. 

As noted above, when making a comparison between the 4FS and 5FS calculations at 
NNLO, we must include in the former the explicit real and virtual 6-quark contributions which 
first appear at 0(a|) in perturbation theory. When calculated with a non-zero 6-quark masqj 
these contributions are finite, and simply add to the 4-fiavour contributions. Sample Feynman 



^The 5FS cross sections are identical to those reported in Ref. [I]. 

^For consistency, we use mi, = 4.75 GeV, the same value used in the 5FS [T] to generate the 6-quark PDFs 
via the DGLAP equation. 
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(a) 



(b) 
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(c) 



(d) 




^Dooom 



^DOOOD^ 



WW 



(e) 



(f) 



Figure 11: Sample Feynman diagrams for the various (a,b,c,d) virtual and (e,f) real 6-quark contri- 
butions to the Z cross section. The 6-quarks are shown as the thicker fermion lines. 
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subprocess 


AfoO"^ (Tevatron) 


Abcx^ (LHC, 7 TeV) 


Aba^ (LHC, 14 TeV) 


q + q^ Z 


5.230 X 10-6 


-2.124 X 10-5 


-6.440 X 10-5 


q + q^ Z + g 


4.901 X 10"^ 


6.185 X 10-5 


9.701 X 10-5 


q{q) +g^Z + q{q) 


-2.862 X 10-5 


-1.456 X 10-4 


-2.632 X 10-4 


q + q^ Z + b + b 


3.754 X 10-4 


1.450 X 10-3 


3.382 X 10-3 


g + g ^ Z + b + b 


2.090 X 10-4 


5.287 X 10-3 


1.997 X 10-2 


total 


6.100 X 10-4 


6.632 X 10-3 


2.312 X 10-2 



Table 8: Additional 0(a;|) contributions to the total Z 4FS NNLO cross section in nb (multiplied 
by leptonic branching ratio) at the Tevatron and LHC arising from real and virtual 6-quark processes. 



diagrams for the various real and virtual 6-quark contributions are shown in Fig. [TTl and 
analytic expressions are given in Refs. [83|l8^ . We summarise these in Appendix \^ and derive 
small-mass (i.e. ml/M^ ^ 1) expansions which are useful in practice. The results are presented 
in Table [HI for Tevatron and {y/s = 7, 14 TeV) LHC energies. Evidently the 2 — )■ 3 process 
contributions, Eq. (IA.7I) . are by far the most dominant, and especially so at the LHC where 
they are two orders of magnitude larger then the rest of the 0(a|) contributions combined. 
In Table [9] we add these real and virtual 6-quark contributions to the bulk of the NNLO cross 
section coming from 4FS light quarks and gluons, and compare the total with the benchmark 
5FS NNLO results presented in Ref. [1], in which the b quark is treated as a massless parton in 
the subprocess cross sections, i.e. both in the initial and final states and in loop contributions. 
Also shown in Table |9] are the {bb, bg, . . .) contributions in which the Z couples directly to a 
b quark and in which there is at least one 6 or 6 quark in the initial stateo These represent 
0.03%, 1.5% and 2.4% of the total 5FS NNLO cross sections at the Tevatron, LHC(7 TeV) and 
LHC (14 TeV) respectively. 

From Table HI we see that at the Tevatron the 5FS and 4FS-|-6-quark total cross sections 
are the same to within 1%. At the LHC, however, the additional 6-quark contributions to the 
4FS cross section are at the +1% level and so do not completely compensate the "missing" 
bb ^ Z contributions to the 5FS cross section. To be precise. 



'NNLO 



(5FS,6) = 48.4 (14.5) pb > A^a^ = 23.1 (6.6) pb 



(8) 



at ^/s = 14 TeV (7 TeV), which results in the "full" 4FS total cross section being 2.3% 
(1.7%) smaller than the the 5FS cross section. We interpret this as meaning that the DGLAP- 
resummed [a5'ln(M|/m^)]"' contributions that are absorbed into the b parton distribution are 
numerically important. 

Our results are in broad agreement with the study of Ref. [85], in which the different 
calculational methods for heavy particles produced in association with b quarks, gg — )■ Xbb 
(4FS) versus bb ^ X (5FS), were studied in detail for X = H, Z. In particular, it was shown 



^"This excludes small 0(0^) contributions initiated by light quarks and gluons, e.g. qq,gg 
the Z couples to b quarks. 



Zbb, in which 
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Tevatron, s/s = 1.96 TeV 


i?-a|NLo(4FS)(nb) 


5-a|NLo(5FS)(nb) 


i?-a|NLo(5FS,6)(nb) 


^0^ 

a! 


0.2013 
0.0409 
0.0063 


0.2016 
0.0431 
0.0060 


0.0012 
-0.0002 
-0.0003 


total 


0.2485 


0.2507 


0.0008 


A,a^ 


0.0006 


— 




total + Afecr^ 


0.2491 


0.2507 






LHC, v^ = 7 TeV 


fi-a|NLo(4FS)(nb) 


i?-4i,Lo(5FS)(nb) 


5.a|NLo(5FS,6)(nb) 


^0^ 


0.8083 
0.1239 
0.0037 


0.8266 
0.1322 
-0.0002 


0.0202 
-0.0020 
-0.0037 


total 


0.9359 


0.9586 


0.0145 


A,a^ 


0.0066 


— 




total + AfeO-^ 


0.9426 


0.9586 






LHC, v^ = 14 TeV 


fi-a^NLo(4FS)(nb) 


i?-a|j,Lo(5FS)(nb) 


i?-4NLo(5FS,6)(nb) 


^0^ 


1.7472 
0.2384 
-0.0047 


1.8110 
0.2557 
-0.0153 


0.0641 
-0.0050 
-0.0107 


total 


1.9809 


2.0514 


0.0484 


A,a^ 


0.0231 


— 




total -f- Aftcr^ 


2.0040 


2.0514 





Table 9: NNLO predictions for the total Z cross section (multiplied by leptonic branching ratio 
B) at the Tevatron and LHC using MSTW 2008 NNLO PDFs [Ij as input, broken down into the 
ag {n = 0,1,2) contributions, with {q = 
the 4FS calculation and {q = u,d,s,c,b; 



LXi t, Ci • o . C-. 



a 



(5). 



af- 4-flavour MSTW 2008 NNLO PDFs} in 
5-flavour MSTW 2008 NNLO PDFs} in the 5FS 
calculation. The final column gives the contribution to the 5FS cross sections from processes where 
the Z couples directly to h quarks. The additional 0{ag) contributions to the cross section arising 
from real and virtual 6-quark processes, taken from Table |H1 are added to the 4FS cross section in 
the last line of each sub-table. 
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that with the canonical choice of scale /i = M^, the LO gg -^ Zbb cross section at the LHC 
was a factor of two smaller than the NNLO bb ^- Z cross section, consistent with our results 
in Eq. ([H]) above. It was also shown in Ref. [HS] that the agreement between the 4FS and 5FS 
calculations improves if the scale fi is reduced to around Mz/3 (found by choosing fi ~ ^/^ 
near the end of the collinear plateau in the quantity —tda/dt for the process gb — t- Zb). The 
NNLO bb ^ Z cross section is approximately scale independent, while the LO gg ^^ Zbb cross 
section increases with decreasing scale, primarily because of the overall q;| factor. 

Of course the explicit gg — )■ Zbb contribution corresponds only to the n = 1 term in the 
resummed [a5ln(M|/m^)]"' perturbation series implicit in the bb ^ Z 5FS calculation. Com- 
plete agreement between the two schemes would be obtained in a fully all-orders perturbative 
calculation. Note that since at all collider energies Af,a^ is dominated by the contributions 
involving the Zbb final state, we would expect that higher-order corrections to the Zbb produc- 
tion process, which is here calculated only at leading order, will generate the [as'ln(M|/m^)]" 
terms implicit in the 6-quark PDF. The NLO (i.e. 0(a|)) corrections to the Zbb total cross 
sections with m^ 7^ have recently been calculated [86H891. Although the results presented in 
Ref. |89] impose a minimum p^ (> 5 GeV), the i^-factoio is evidently rather independent of 
p\, at small p^, suggesting that the NLO/LO i^-factor for the fully inclusive Zbb cross section is 
approximately 1.5 for the LHC at y/s = 14 TeV. It is therefore plausible that even higher-order 
perturbative corrections can account for the factor of 2 difference in the 4FS and 5FS cross 
sections, Eq. ([8]). This conclusion is supported by the fact that the scale dependence of the 
4FS calculation for Zbb production at NLO is only mildly weaker than at LO [861489] . 

In Ref. |9Tj it was shown that for x in the region of 0.01-0.05, the relevant region for Z + bb 
production at the 14 TeV LHC, the ratio of the GM-VFNS structure function F^ to the FENS 
structure function was ~ 1.5 at LO at high scales. This represents the effect of either resumming 
the [as ln(Q^/'m^)]" contributions, or keeping only the contribution of fixed order in as for one 
parton. For hadron-hadron processes we would expect the difference to be about 1.5^ > 2, 
exactly as observed. At NLO for structure functions at this x the ratio is reduced to ~ 1.1, 
so inclusion of the extra ln((5^/m^) removes much of the discrepancy present at LO. However, 
for hadron-hadron processes, NLO in the fixed flavour scheme only contains the extra large 
logarithm for one of the two incoming partons, so the ratio between the 5FS and 4FS would be 
roughly 1.5 x 1.1 ^ 1.6, again as we expect to see in practice. It would only be at NNLO in the 
4FS, when the double logarithm for both incoming PDEs is included, that we would expect to 
see the reduction to roughly 1.1^ ^ 1.2 in the ratio of the 5FS to 4FS cross sections. This is a 
general feature, i.e. the 4FS (or 3FS) will converge to the resummed 5FS results more slowly 
for hadronic processes than for those in DIS. 

In summary, the 5FS PDEs are clearly the most appropriate to use for inclusive quantities 
such as a^ (or cr'^'J''*, etc.) at high Q^ scales, where resummation of the [ag ln((5^/'n2^)]"' 
contributions is evidently important. However, for more exclusive quantities like a^^'', where 



-'^^The iiT-factor is here defined as the ratio of cross sections calculated using the dynamical scale fi^ 
Ml + (p^^)2 + (p^^)2, and with the standard 5-fiavour CTEQ6M/CTEQ6L1 PDFs W at NLO/LO. 
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the b quarks are measured in the final state, the 4FS parton distributions are more appropriate 
since the 5FS calculation gives no information on the spectator b and b quarks which must be 
present in the final state for the b- and 6-initiated processes. Note that if only the total cross 
section is required, without cuts imposed on the 6-quarks, then a 5FS is still better, e.g. for Zbb 
a 5FS calculation can be used for 66 — )■ Z at NNLO, where the 6-quarks couple directly to the 
Z, and so there are implicitly also two 6-quarks in the final state [BS]- However, if cuts must be 
applied to the 6-quarks, as is the case in the experimental measurement, then a 4FS calculation 
is more appropriate. Similar remarks apply to the calculation of Hbb production [92149^ and 
other processes where 6-quarks are detected in the final state. In a recent study [95], the 
production of a charged Higgs boson in association with a t quark was considered in both the 
4FS {gg — )■ tbH~ etc.) and 5FS {gb — )■ tH~ etc.) to NLO in pQCD, using the appropriate 
MRST 2004 PDF sets [OlITU] . The central predictions in the 5FS were shown to be approximately 
40% larger than those in the 4FS. Even taking the scale uncertainty into account the 4FS and 
5FS NLO cross sections are barely consistent. 

An ideal calculation would combine the best features of the 4FS and 5FS so as to resum 
[aslTo.{Q'^/7nl)]^ terms while also retaining a finite rrib dependence in the partonic cross section 
(rather than setting mj, to zero as done in the 5FS). This matching has, of course, been done 
for structure functions in DIS using different variants of the GM-VFNS, but applications to 
hadron collider cross sections are more involved and have so far been limited (see Ref. [nS] for an 
application of the GM-VFNS to the pt spectrum in heavy-fiavour hadroproduction). However, 
for processes where the hard scale is, for example, Q^ ~ M|, then the GM-VFNS calculation 
will differ from the ZM-VFNS (5FS) only by terms 0{ml/Ml ~ 0.3%). We would therefore 
expect the complete GM-VFNS calculation to give results very close to the pure ZM-VFNS 
(5FS) for the total cross section. 

Note that rather than producing separate 4-flavour PDFs for use in a 4FS calculation, an 
alternative approach (e.g. used at NLO in Refs. [971[98]) is to use the conventional 5-fiavour 
PDFs, then pass to the 4-fiavour scheme using counterterms given in Ref. |96] . However, these 
counterterms [96] are equivalent to using the inverse of transition matrix elements, but only 
out to order as- One could indeed use the transition matrix elements themselves to go from 
4-fiavour to 5-fiavour PDFs, except this would not sum the logarithmic terms, [as ln((5^/m^)]", 
in the PDF evolutiono Hence, the use of counterterms [96] is a less complete way of going 
from a 5FS to a 4FS, and instead, we recommend that dedicated 4-fiavour PDFs be used in 4FS 
calculation's Previously, a major advantage ot using 5-fiavour PDFs with counterterms was 
that eigenvector PDF sets to calculate PDF uncertainties were not made available for existing 



^^See Fig. 9 of Ref. _80; for a comparison of 5-flavour NNLO PDFs obtained from 3-flavour PDFs either by 
evolution or by applying fixed-order matching conditions; the differences will be larger at NLO. 

^^However, for the 4FS calculation of i-channel single-top production at the Tevatron and LHC |97|, it was 
explicitly checked that results obtained with the dedicated 4-flavour MRST set [9] were consistent (within 
the numerical integration precision) with those obtained with the corresponding 5-flavour MRST set [10] plus 
appropriate counterterms !96J. We thank R. Frederix and F. Tramontano for discussions on this issue. 
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4-flavour PDFs [9] . However, we have now provided eigenvector PDF sets also for the 4-flavour 
PDFs, therefore this advantage no longer holds. 



5 Conclusions 

We have repeated the NLO and NNLO MSTW 2008 global PDF analyses [Ij for a range of 
heavy-quark masses about their "default" values rric = 1.40 GeV and nib = 4.75 GeV. For the 
charm quark, we found that the global data prefer the values rric = 1.45 (1.26) GeV at NLO 
(NNLO). The most discriminating data are, as anticipated, the HERA data for F2 [HHSO] . On 
the other hand, for the bottom quark, the data included in the global fit (excluding ^2^) do 
not put a meaningful constraint on the value of rrib, while the HERA F2 data slightly favour 
rub ~ 4.75-5 GeV. We pointed out that precise determinations of the heavy-quark masses in 
the MS scheme are affected by poorly convergent perturbative series in the conversion to the 
pole masses, particularly for the case of the charm quark. Recent precise combined HERA data 
on a^'^ [79] and F2 [71] will in future be able to narrow the favoured range of the charm-quark 
pole mass rric- Note, however, that uncertainties from the choice of GM-VFNS [61] mean that 
the favoured value of rric will be correlated to some extent with the particular choice of GM- 
VFNS, although this correlation will be much smaller at NNLO than at NLO, as will other 
uncertainties arising from the choice of GM-VFNS |64j . 

We explored the effect of the values of the heavy-quark masses on W, Z and Higgs pro- 
duction at the Tevatron and LHC. Varying the charm mass by ±0.15 GeV changes the cross 
sections by ±1% or less at Tevatron energies and by ±2% at the LHC energy of y/s = 14 TeV. 
The various weak boson cross-section ratios are essentially unchanged. The predictions for W, 
Z and Higgs cross sections are much less dependent on the value taken for m^. We provided a 
recommendation on how to include the uncertainty arising from the choices of rric and rrib in a 
generic cross-section calculation. 

We also presented PDF sets obtained in a framework with different active numbers of quark 
flavours, as done previously in the context of the MRST 2004 analysis [H]. Explicitly, we 
determined 4-fiavour PDF sets in which charm becomes an active parton, but bottom does not, 
and 3-fiavour PDF sets where charm and bottom are not partons, but only appear in the final 
state. The analogous 5-flavour parton sets are simply those of MSTW 2008 [Ij. Of course, the 
latter, which in the absence of top corresponds to PDFs of a (general-mass) variable flavour 
number scheme, are generally to be preferred, particularly for inclusive quantities at high Q^ 
scales where the resummation of the [asin.{Q'^/ml)]'^ contributions is essential. However, for 
more exclusive processes, such as Zbb, Hbb, . . ., where 6-quarks are observed in the final state, 
the 4-fiavour parton distributions are more appropriate. For illustration, we computed the 
various components of the Z production cross section to 0(«|) at the Tevatron and LHC, and 
compared the predictions obtained using the 4FS and 5FS (= MSTW 2008 [Ij) parton sets. 
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The additional grids for all PDF sets discussed in this paper are made publicly available [60] , 
for use either with the standalone MSTW interpolation code or via the lhapdf interface PU] . 
To be precise, grids for the following PDF sets are made available: 

• For the default quark masses, rric = 1.40 GeV and rrib = 4.75 GeV, we provide LO, NLO 
and NNLO grids for 3- and 4-flavour PDFs (central set and 40 eigenvector sets at both 
68% and 90% C.L.). These grids complement the existing grids for the 5-flavour PDFs. 

• For rric in the range 1.05 GeV to 1.75 GeV (in steps of 0.05 GeV), we provide NLO and 
NNLO grids for 3- and 5-flavour PDFs (central set only) for both free and fixed as{M^). 

• For rrih in the range 4.00 GeV to 5.50 GeV (in steps of 0.25 GeV), we provide NLO and 
NNLO grids for 4- and 5-fiavour PDFs (central set only) for free as{M'^) only. 

These additional grids should prove to be useful in future for detailed studies of a variety of 
collider processes involving heavy quarks. 



A Appendix 

We consider the 6-quark contributions to the 4FS calculation, closely following the discussion 
in Ref. [HS]. First, we have the two-loop corrections to the g + g — )■ Z Born process, examples 
of which are shown in Fig. [TTlfa) and (b). In the notation of Ref. [83] . 



W^, = 6(1 - r)C,T, (^) 



2agaQGi{pb) + -{vq+ aj) F{pb 



(A.l) 



with pb = ml/Mz ~ 2.7 x 10 ^. The functions F and Gi are defined in Refs. [83] and [S 
respectively. Since in practice p^, ^ 1, we can use the following small-mass expansions: 

^, , 4 „ 38,2 A6, 530\ , 3355 152, 16, 

^'^^ 9 9 V 3 27/81 9 ^ 3 ^ 

+p i 16L^ - 32L + 128 - —(2 ] + C(p'), (A.2) 

Gi(p) ^ 4p(-2C2L + 2C3-3) + 0(p2), (A.3) 

with L = ln(l/p) ^ 5.9 for p = pb- The contributions to the total Z cross section from these 
loop diagrams are given in the first row of Table [HI They are numerically very small and 
negative (positive) at the LHC (Tevatron)o 



^*Even though the function G\{p) is suppressed by a power of p relative to F{p) at small p, it gives a 
comparable contribution to the cross section. This is because F has a zero close to the physical value p — pb- 
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Next we have the one-loop corrections to the 2—7-2 processes q + q ^^ Z + g and g + q{q) — )■ 
Z + g(g), examples of which are shown in Fig. ITTT c) and (d): 



_ 1 

qq- 2' 



W^f- = ^raqaqCpTf 



as_ 

IT 



1 + f 
1-f 



{-2 + 2f ( Ji(4p5f) - Ji(4pfe))} - 4pbr (J2(4p6f) - J2{^Pb)) 



K = l<^q<^QTf{^yH,if, 



Pb) 



(A. 
(A.5) 



where s = f ^M| and the functions Hi, Ji and J2 are defined in Ref. |84j. The following 
small-mass approximations are again useful: 



Hi{t,0) = 2r 

Mp) ' 

J2{p) 



2{t - 1) + In 



P 



In 



1-r 



r 



Li2 



r-1 



r 



2(1 -r), 



L + ln4 - 2 - ^ (L + ln4 + 1) + 0(p2), 



P 



-((L + ln4f-7r^)-|(L + ln4) + 0(/). 



Finally, we have the 2—7-3 processes 



q + q^ Z + h + h, g + g^Z + b + b, 



(A.6) 



(A.7) 



see Fig. [TlT e) and (f). The calculation of these is in principle straightforward, with matrix 
elements squared integrated over three-body phase space. The non-zero 6-quark mass makes 
these contributions infra-red and collinear finite, i.e. the singularities present for massless quarks 
are replaced by terms with logarithmic (ln(M|/m^)) behaviour in the small rrib limit. We use 
the MCFM implementation [100] of the matrix elements and phase space, with exactly the same 
parameter choice as in the previous MSTW calculations. The results are given in Table [H Note 
that the gg contribution is larger than the qq contribution at the LHC, while the converse is 
true at the Tevatron. 



References 

[1] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 63 (2009) 189 
[arXiv:0901.0002l [hep-ph]]. 

[2] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Phys. Lett. B 652 (2007) 292 
[a rXiv:0706.0459l [hep-ph]]. 

[3] H. L. Lai, M. Guzzi, J. Huston, Z. Li, P. M. Nadolsky, J. Pumplin and C. P. Yuan, Phys. 
Rev. D 82 (2010) 074024 [ar Xiv: 1007.22411 [hep-ph]]. 

[4] P. M. Nadolsky et al, Phys. Rev. D 78 (2008) 013004 |arXiv:0802.0007l [hep-ph]]. 



34 



[5] R. D. Ball, L. Del Debbio, S. Forte, A. Guffanti, J. I. Latorre, J. Rojo and M. Ubiali, 
Nucl. Phys. B 838 (2010) 136 |arXiv: 1002.44071 [hep-ph]]. 

[6] A. D. Martin, W. J. Stirling, R. S. Thorne and G. Watt, Eur. Phys. J. C 64 (2009) 653 
[arXiv:0905.353T1 [hep-ph]]. 

[7] H. L. Lai, J. Huston, Z. Li, P. Nadolsky, J. Pumplin, D. Stump and C P. Yuan, Phys. 
Rev. D 82 (2010) 054021 [i ^XTv: 1004.4624 [hep-ph]]. 

F. Demartin, S. Forte, E. Mariani, J. Rojo and A. Vicini, Phys. Rev. D 82 (2010) 014002 
[arXiv: 1004.09621 [hep-ph]]. 

[9] A. D. Martin, W. J. Stirling and R. S. Thorne, Phys. Lett. B 636 (2006) 259 
[arXiv:hep-ph/0603143 |. 

[10] A. D. Martin, R. G. Roberts, W. J. Stirhng and R. S. Thorne, Phys. Lett. B 604 (2004) 
61 [ar Xiv:hep-ph/0410230] . 

[11] M. Buza, Y. Matiounine, J. Smith and W. L. van Neerven, Eur. Phys. J. C 1 (1998) 301 
|arXiv:hep-ph /9612398]. 

[12] I. Bierenbaum, J. Blumlein and S. Klein, Nucl. Phys. B 820 (2009) 417 |arXiv:0904.3"563l 
[hep-ph]]. 

[13] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B 392, 162 
(1993). 



[14] R. S. Thorne, Phys. Rev. D 73 (2006) 054019 |arXiv:hep-ph/0601245 



[15] S. Catani, M. Ciafaloni and F. Hautmann, Nucl. Phys. B 366 (1991) 135. 

[16] E. Laenen and S. O. Moch, Phys. Rev. D 59 (1999) 034027 |arXiv:hep-ph/9809550| . 

[17] S. Alekhin and S. Moch, Phys. Lett. B 672 (2009) 166 |arXiv:0811.i4T2] [hep-ph]]. 

[18] N. A. L. Presti, H. Kawamura, S. Moch and A. Vogt, PoS DIS2010 (2010) 163 
[arXiv: 1008.0951. [hep-ph]]. 

[19] B. W. Harris and J. Smith, Nucl. Phys. B 452 (1995) 109 |arXiv:hep-ph/9503484| ; 



[20] B. W. Harris and J. Smith, Phys. Rev. D 57 (1998) 2806 |arXiv:hep-ph/9706334] . 
[21] S. Frixione, P. Nason and B. R. Webber, JHEP 0308 (2003) 007 |arXiv:hep-ph /0305252]. 
[22] A. C. Benvenuti et al. [BCDMS Collaboration], Phys. Lett. B 223 (1989) 485. 
[23] A. C. Benvenuti et al. [BCDMS Collaboration], Phys. Lett. B 237 (1990) 592. 

35 



[24] M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 483 (1997) 3 



|arXiv:hep-p h/961023ll . 

[25] M. Arneodo et al. [New Muon Collaboration], Nucl. Phys. B 487 (1997) 3 



[arXiv:hep-ex/9611022| . 

[26] M. R. Adams et al. [E665 Collaboration], Phys. Rev. D 54 (1996) 3006. 

[27] L. W. Whitlow, E. M. Riordan, S. Dasu, S. Rock and A. Bodek, Phys. Lett. B 282 (1992) 
475. 

[28] L. W. Whitlow, Ph.D. thesis, Stanford University, 1990, SLAC-0357. 

[29] L. W. Whitlow, S. Rock, A. Bodek, E. M. Riordan and S. Dasu, Phys. Lett. B 250 (1990) 
193. 



[30] J. C. Webb, Ph.D. thesis. New Mexico State University, 2002, |arXiv:hep-ex/030103l] 



Paul E. Reimer, private communication (for the radiative corrections). 

[31] R. S. Towell et al. [FNAL E866/NuSea Collaboration], Phys. Rev. D 64 (2001) 052002 
[arXiv:hep-ex/0103030i . 

[32] M. Goncharov et al. [NuTeV Collaboration], Phys. Rev. D 64 (2001) 112006 



(arXiv:hep-ex/0102049) . 

[33] M. Tzanov et al. [NuTeV Collaboration], Phys. Rev. D 74 (2006) 012008 
[a rXiv:hep-ex/0509010] . 

[34] G. Onengut et al. [CHORUS Collaboration], Phys. Lett. B 632 (2006) 65. 

[35] E. M. Lobodzinska [HI Collaboration], larXiv:hep-ph/0311180l 



[36] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 21 (2001) 33 [arXiv:hep-ex/0012053|. 



[37] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 19 (2001) 269 [arXiv:hep-ex/0012052|. 



[38] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 30 (2003) 1 [ arXiv:hep-ex/0304003| . 

[39] J. Breitweg et al. [ZEUS Collaboration], Eur. Phys. J. C 7 (1999) 609 
|arXiv:hep-ex /9809005] . 

[40] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 21 (2001) 443 
[a rXiv:hep-ex/0 105090| . 

[41] S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 28 (2003) 175 
[arXiv:hep-ex/0208040j . 



36 



[42] S. Chekanov et al. [ZEUS Collaboration], Pliys. Rev. D 70 (2004) 052001 
[arXiv:hep-ex/0401003| . 

[43] S. Chekanov et al. [ZEUS Collaboration], Eur. Pliys. J. C 32 (2003) 1 
[a rXiv:hep-ex/0307043] . 

[44] C. Adloff et al. [HI Collaboration], Z. Pliys. C 72 (1996) 593 ||arXiv:hep-ex/9607012| . 



[45] C. Adloff et al. [HI Collaboration], Pliys. Lett. B 528 (2002) 199 |arXiv:hep-e?^1 08039 ] . 

[46] A. Aktas et al. [HI Collaboration], Eur. Pliys. J. C 45 (2006) 23 [ arXiv:hep-ex/0507081| . 

[47] A. Aktas et al. [HI Collaboration], Eur. Pliys. J. C 40 (2005) 349 |arXiv:hep-ex /0411046|. 

[48] J. Breitweg et al. [ZEUS Collaboration], Eur. Pliys. J. C 12 (2000) 35 
|a rXiv:hep-ex/9908012] . 

[49] S. Chekanov et al. [ZEUS Collaboration], Phys. Rev. D 69 (2004) 012004 
(a^v:hep-ex/0308068| . 



[50] S. Chekanov et al. [ZEUS Collaboration], JHEP 0707 (2007) 074 [arXiv:0704.3"562l [hep- 

ex]]. 

[51] A. Aktas et al. [HI Collaboration], Phys. Lett. B 653 (2007) 134 [ar2£iviQ70O722] [hep- 
ex]]. 

[52] S. Chekanov et al. [ZEUS Collaboration], Phys. Lett. B 547 (2002) 164 
(a^v:hep-ex/0208037| . 

[53] S. Chekanov et al. [ZEUS Collaboration], Nucl. Phys. B 765 (2007) 1 
[arXiv:hep-ex/0608048 1 . 

[54] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. Lett. 101 (2008) 062001 
|arXiv:0802.2"400l [hep-ex]]. 

[55] A. Abulencia et al. [CDF - Run II Collaboration], Phys. Rev. D 75 (2007) 092006 
[Erratum-ibid. D 75 (2007) 119901] [arXiv:hep-ex/0701051j . 

[56] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. D 77 (2008) 011106 [arXiv:0709.4"254l 

[hep-ex]]. 

[57] D. Acosta et al. [CDF Collaboration], Phys. Rev. D 71 (2005) 051104 
[aiOGvrhep-ex/0501023] . 



[58] V. M. Abazov et al. [D0 Collaboration], Phys. Rev. D 76 (2007) 012003 
[arXiv:hep-ex/0702025j . 



37 



[59] J. Han et al. [CDF Collaboration], "dcr/dy distribution of Drell-Yan dielectron pairs", 



[6o; 

[61 

[62 

[63 

[64 
[65 
[66 

[67; 

[68 

[69 

[70 
[71 

[72 
[73 
[74; 

[75 
[76 
[77 



Public Note, May 2008, http : //www-cdf . f nal . gov /phy sics/ewk/2008/dszdy/ 
http : //pro j ects . hepf orge . org/mstwpdfT] 



F. D. Aaron et al. [HI Collaboration], Eur. Phys. J. C 65 (2010) 89 |arXiv:0907.2643l 

[hep-ex]]. 

H. Abramowicz et al. [ZEUS collaboration], Eur. Phys. J. C 69 (2010) 347 
|arXiv:1005.3 396 [hep-ex]]. 

S. Chekanov et al. [ZEUS Collaboration], Eur. Phys. J. C 65 (2010) 65 |arXiv:0904.3487l 

[hep-ex]]. 



R. S. Thorne, PoS DIS2010 (2010) 053 [a rXiv:1006.59251 [hep-ph]]. 

C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1. 

K. G. Chetyrkin, J. H. Kuhn, A. Maier, P. Maierhofer, P. Marquard, M. Steinhauser and 
C. Sturm, Phys. Rev. D 80 (2009) 074010 [arXiv:0907.2TT0l [hep-ph]]. 

S. Narison, Phys. Lett. B 693 (2010) 559 [ arXiv: 1004. 53331 [hep-ph]]. 

K. G. Chetyrkin and M. Steinhauser, Nucl. Phys. B 573 (2000) 617 
arXiv:hep-ph/9911434| . 



K. Melnikov and T. v. Ritbergen, Phys. Lett. B 482 (2000) 99 |arXiv:hep-ph/9912391 

J. H. Kuhn, Nucl. Phys. Proc. Suppl. 181-182 (2008) 141. 

I. I. Y. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vainshtein, Phys. Rev. D 50 (1994) 
2234 ||arXiv:hep-ph/9402360|. 



M. Beneke and V. M. Braun, Nucl. Phys. B 426 (1994) 301 [arXiv:hep-ph/9402364| . 



A. H. Hoang and A. V. Manohar, Phys. Lett. B 633 (2006) 526 |arXiv:hep-ph /0509195]. 

HERA Heavy Flavour Combination Group, "Combination of Fg'^ from DIS measurements 
at HERA", Hlprelim-09-171, ZEUS-prel-09-015. 

J. J. Aubert et al. [European Muon Collaboration], Nucl. Phys. B 213 (1983) 31. 

A. M. Cooper-Sarkar, PoS DIS2010 (2010) 023 |arXiv:1006.4"47n [hep-ph]]. 

M. Gluck, P. Jimenez-Delgado and E. Reya, Eur. Phys. J. C 53 (2008) 355 
[arXiv:0709.0614 [hep-ph]]; 



[78] P. Jimenez-Delgado and E. Reya, Phys. Rev. D 79 (2009) 074023 |arXiv:0810.4274l [hep- 
ph]]. 



38 



[79] F. D. Aaron et al. [HI and ZEUS Collaborations], JHEP 1001 (2010) 109 
|arXiv:0911.0 884 [hep-ex]]. 

[80] S. Alekhin, J. Blumlein, S. Klein and S. Moch, Pliys. Rev. D 81 (2010) 014032 
|arXiv:0908.2766l [hep-ph]] . 

^1] R. S. Thorne and W. K. Tung. ^ arXiv:0809.0714l [hep-ph]. 

[82] M. Buza and W. L. van Neerven, Nucl. Phys. B 500 (1997) 301 |arXiv:hep-ph/9702242| . 

[83] P. J. Rijken and W. L. van Neerven, Phys. Rev. D 52, 149 (1995) | arXiv:hep-ph/9501373| . 

[84] R. J. Gonsalves, C. M. Hung and J. Pawlowski, Phys. Rev. D 46 (1992) 4930. 

[85] F. Maltoni, T. McElmurry and S. Willenbrock, Phys. Rev. D 72 (2005) 074024 
[arXiv:hep-ph/0505014] . 

[86] F. Febres Cordero, L. Reina and D. Wackeroth, Phys. Rev. D 74 (2006) 034007 
[arXiv:hep-ph/0606102] ; 

F. Febres Cordero, L. Reina and D. Wackeroth, Phys. Rev. D 78 (2008) 074014 
|arXiv:0806.0 808 [hep-ph]]; 

F. Febres Cordero, L. Reina and D. Wackeroth, Phys. Rev. D 80 (2009) 034015 
[a rXiv:0906.1923] [hep-ph]]. 

[89] F. Febres Cordero and L. Reina, "Z66 for CMS, Study I", private communication. 

[90] J. Pumphn, D. R. Stump, J. Huston, H. L. Lai, P. M. Nadolsky and W. K. Tung, JHEP 
0207 (2002) 012 |arXiv:hep-ph/0201195| . 



[91] R. S. Thorne and R. G. Roberts, Phys. Rev. D 57 (1998) 6871 |arXiv:hep-ph/9709442|. 



[92] S. Dittmaier, M. Kramer and M. Spira, Phys. Rev. D 70 (2004) 074010 
|arXiv:hep-ph/0309204] . 



[93] S. Dawson, C. B. Jackson, L. Reina and D. Wackeroth, Phys. Rev. D 69 (2004) 074027 
|arXiv:hep-ph/03110'67] . 



[94] J. M. Campbell et a/., arXiv:hep-ph/0405302 



[95] S. Dittmaier, M. Kramer, M. Spira and M. Walser. [arXiv:0906.2648l [hep-ph]. 

[96] M. Cacciari, M. Greco and P. Nason, JHEP 9805 (1998) 007 |arXiv:hep-ph79803400] . 

[97] J. M. Campbell, R. Frederix, F. Mahoni and F. Tramontano, Phys. Rev. Lett. 102 (2009) 
182003 ^arXiv:0903.0 005 [hep-ph]]. 



39 



[98] J. M. Campbell, R. Frederix, F. Maltoni and F. Tramontano, JHEP 0910 (2009) 042 
|arXiv:0907.3 933 [hep-ph]]. 



[99] M. R. Whalley, D. Bourilkov and R. C. Group, arXiv:hep-ph/0508110, 



LHAPDF code from http://projects.hepforge.org/lhapdf/ 



[100] J. M. Campbell and R. K. Ellis, Phys. Rev. D 62 (2000) 114012 [a rXiv:hep-ph/0006304| ; 
MCFM code from http : //mcf m . f nal . gov/ 1 



40 



